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The morphology of plastic wires packed into spherical confinements is investigated. The total 
length of the injected wire exhibits a power-law scaling with the system size (the ratio between 
container and wire radii R/r), with an exponent dependent on the wire-wall friction. At high friction 
regime, coils rarely form and the number of folds follows a power-law dependence on R/r. During 
the injection process, the fold-size distribution gradnally broadens and becomes more asymmetric, 
reflecting the increase of spatial exclusion effects. The cumulative wire length vs the fold size 
collapses onto a universal curve, allowing for the prediction of the total length of crumpled wire for 
different R/r, material properties, and insertion forces. A self-avoiding random walk model for such 
a far-from-equilibrium process is introduced which remarkably reproduces the experimental data. 

PACS numbers: 05.40.-a, 62.20.F-, 89.75.Da, 46.65.-|-g 


Crumpling of thin objects in confined geometries is 
ubiquitously observed in nature. Examples range from 
DNA packing in viral and bacteriophage capsids Qi 
to folding of insect wings in cocoons Q and plant leaves 
in buds [^. In response to the applied forces and con¬ 
straints, complex morphologies may be adopted where 
friction, nonlinear response, self-avoiding interactions, fi¬ 
nite size effects, and boundary conditions play crucial 
roles. There have been numerical efforts to identify the 
separate contribution of self-avoidance on the morpho¬ 
logical properties of crumpled objects It was found 
that the exponent of the power-law governing the mass- 
size relation differs between phantom elastic sheets and 
when self-avoidance and plasticity are considered d, Q . 
Another numerical study revealed that self-avoidance in 
highly compressed elastic papers leads to self-correlations 
which modify the hierarchical folding process Q. It 
is, however, difficult to disentangle the effects of self¬ 
avoidance and other key factors on the resulting struc¬ 
tures of frictional elasto-plastic materials in experiments. 


So far, very little is known about the structural evo¬ 
lution during crumpling isa. Understanding how the 
space-filling properties of crumpled objects develop is of 
particular importance in technological applications such 
as endovascular coiling treatment of cerebral aneurysms 
[ 1 ^ . As expected from far-from-equilibrium processes, 
the imposed constraints and initial conditions cannot 
uniquely determine the crumpled state. There is indeed 
an ensemble of admissible configurations, from which 
some structural properties of the system can be derived. 


In this Letter, we study the injection of solder wires 
into frictional rigid confinements, which enables us to iso¬ 
late the influence of self-avoiding interactions. By mon¬ 
itoring the morphological evolution, a gradual crossover 
from random to correlated folding events is observed due 
to increasing spatial exclusion effects. The folding pro¬ 


cess of plastic wire can be considered as a self-avoiding 
random walk (SAW) in a confined geometry. While nu¬ 
merical sampling of the ensemble of such SAW trajec¬ 
tories is unfeasible via stochastic Markovian dynamics 
approaches, employed e.g. in the simple case of regular 
lattices, we introduce an alternative approach which ac¬ 
counts for the time evolution of the fold length. The re¬ 
sults of Monte Carlo simulations show satisfactory agree¬ 
ment with experiments. We verify that the final length 
of the injected wire can be estimated from the geometry, 
material properties, and imposed constraints. 

Setup — The experimental setup consists of a rigid 
hollow spherical container of inner radius R with a small 



FIG. 1. (a) Morphologies of plastic wires packed into spheri¬ 
cal containers, at low (top) and high (bottom) wire-container 
friction (b) The total length L of the inserted wire, 

scaled by the cross-section radius r of the wire versus the 
relative system size R/r. The dashed lines indicate the best 
power-law fits. Inset: Schematic of the side view of the setup. 
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hole to insert the plastic wire (see Fig. [T|). Several poly¬ 
meric semi-transparent molds with radii R ranging be¬ 
tween 5 and 20 mm are constructed. A small nozzle 
and two counterrotating rollers are attached to the in¬ 
jecting hole, which facilitate the control of the insertion 
force and speed. The friction coefficient between the 
wire and the container wall is 0.4 (0.2) when the 

smooth inner wall of each mold is roughened by sand¬ 
paper (lubricated with silicone oil). Traditional solder 
wire Sn 6 oPb 4 o with Young’s modulus Y of around 30 GPa 
and dry wire-wire friction coefficient /i^ ^~0.2 is used as 
working material. We insert wires with cross-section ra¬ 
dius r = 0.4, 0.5, 0.6, or 0.8 mm into the ball through the 
nozzle with feeding speed of about 1 mm/s. The insertion 
process continues until the applied feeding force is unable 
to fold the wire anymore, and the wire tends to buckle 
outside of the ball. In the absence of elastic restoring en¬ 
ergy, the final structure preserves its shape after opening 
the mold, which allows for a detailed structural analysis. 

Friction-dependent structures — We observe distinct 
morphological phases, depending on the wire-wall fric¬ 
tion. At low friction (/i^^,j~0.2), the wire can slide on the 
container surface smoothly, which increases the chance 
of coil formation during the early stages of injection pro¬ 
cess [see Fig. (Ha)]. Even highly ordered coils form at the 
outer layer of the crumpled structure in relatively small 
spheres, which are reminiscent of the structures observed 
in packings of elastic wires 13 However, at higher 


friction (/r^^|j~0.4), the wire randomly bends due to local 
constraints instead of relaxing towards a global minimum 
energy. This effect, together with excluded volume inter¬ 
actions, promote the formation of folds rather than coils. 
By introducing the frequency of sharp changes of curva¬ 
ture as the order parameter, it is recently shown that a 
phase transition to disordered structures occurs as fric¬ 
tion is increased m. We roughened the surface of ball 
to ensure that is far above the transition point, i.e. 
mixed fold-coil structures are avoided and folding is the 
dominant process. The wire-wire friction prevents coil¬ 
ing in the core of the structure. This enables us to focus 
on the growing role of self-avoidance by analyzing the 
folding statistics. 

The influence of friction on crumpling morphology is 
also visible when looking at the total length L of the in¬ 
serted wire. Figure [ijb) shows Ljr versus the effective 
system size R/r. The data-collapse on a common curve 
for different combinations of R and r verifies that R/r 
is an appropriate control parameter to describe the be¬ 
havior. For a given value of R/r, L is expectedly smaller 
for higher since randomly oriented folds effectively 
fill the space with a lower packing fraction than a mixed 
structure of folds and coils at low At both wall 

friction coefficients, the system-size dependence of L/r is 
overally captured well by the power-law 


with a friction-dependent exponent D. We obtain D = 
2.48 (2.62) at high (low) by fitting the experimen¬ 
tal data to Eq. [TJ Exponent of D = 2.5 was previously 
obtained for ordered packing of elastic wires in spheri¬ 
cal cavities and similar power-law scaling (with T)~2.75) 
was reported for packing of plastic wires 17|. Balancing 
the confining energy and self-avoidance in similar sys¬ 
tems gives a mean-field exponent D=2.5 , because 

the lowest (harmonic) approximation of the confining en¬ 
ergy yields an energy density of the order of {R/r)~^, 
while the self-avoidance energy density is proportional to 
(L/r)^(i?/r)“®, as it originates from binary contacts be¬ 
tween wires whose density nearly grows as (fP' where 4> is 
the packing fraction. Analogously, an exponent 2.5 can 
be obtained between the area of an uncrumpled sheet and 
its radius of gyration after crumpling [B| • This suggests 
the independence of the fractal dimension from the way 
the plastic wire is packed into the container. For exam¬ 
ple, one can first construct a sheet by regularly folding 
the plastic wire in a plane and then crumple it. Another 
point is that from the scaling of L it follows that the 
packing fraction (j)=Trr'^L/{^ ttR^) scales as ^ {R/r)^~^ 
with a negative exponent, indicating the decay of density 
with increasing the effective system size. 

Fold-size statistics — We now turn to the analysis of 
folds statistics at high friction regime. The points of 
folding are distinguishable by sharp changes of wire ori¬ 
entation. A minimum threshold of ninety degrees is im¬ 
posed for the turning angle of the wire to be identified 
as a folding point. We cut the wire at the folding points, 
straighten the segments, and measure their length using 
digital calipers. Straightening of the curved segments 
rarely allows for fold lengths £ longer than the container 
diameter, but we checked that £max remains smaller than 
ttR in the absence of coils. 

A key observation is the scaling of the total number of 
folds N with the effective system size R/r. As shown in 
Fig. [2J a power-law relation of the form 

N^{R/rf (2) 

holds with /3~1.86. Notably, the number of steps in a 



R/t 


R/r 


R/r 


L/r ^ {R/r) 


D 


( 1 ) 


FIG. 2. Scaling of (a) the number of folds N, (b) the di¬ 
mensionless mean fold size {£)/R, and (c) the standard de¬ 
viation w/R versus the effective system size R/r. The full 
(open) symbols indicate experimental (simulation) results. 
The dashed lines are power-law fits to the experimental data. 
The experiment is repeated 3 times for each value of R/r. 
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FIG. 3. (color online) (a) Probability distribution of the nor¬ 
malized fold size tjR. Symbols (lines) are obtained from ex¬ 
periments (simulations), (b) Comparison between the tail of 
the data shown in part (a) and the corresponding best overall 
fits to Gaussian, log-normal, and gamma functions. The solid 
lines show the simulation results. 


self-avoiding random walk on a cubic lattice scales in 
a similar way with the size of the visited region, even 
though with a somewhat smaller exponent /3~5/3 [Uj 
(for an ordinary random walk /3 equals 2). Similar scal¬ 
ing laws also hold for the number of coils in 3D ordered 
elastic packing (/3~1.5) [1^, or the number of loops in 
2D pac king of crumpled wire at high plasticity regime 
[l^ fiil . [ 2 J. From the scaling of L and iV with Rjr one 
finds that the normalized mean fold size {() jR follows 
(i?/r)^“^“^-“°-^. The best fit to the standard devia¬ 
tion of £ versus the system size also follows a power-law: 
w/i?~ (i?/r)° ®. The same exponent is reported for the 
roughness of crumpled surfaces 0 which might indicate 
a universality class for roughness exponent in crumpled 
systems independent of dimensionality. 

The probability distribution P{i/R) of the normalized 
fold size is plotted for different values of R/r in Fig. |3l 
P{£/R) becomes more asymmetric with the peak shifting 
towards smaller fold sizes by increasing R/r. At earlier 
stages of the injection process in a given sphere size, the 
injected wire proceeds in the container without interact¬ 
ing with the accumulated wire. Assuming that the wire 
bends at a random point between the injecting hole and a 
contact point at the container surface, the resulting fold 
length is a random variable, symmetrically distributed 
between 0 and £inax- By increasing the length of accumu¬ 
lated wire, spatial exclusion effects grow and the injected 
wire cannot easily proceed through the sphere without 
touching the crumpled structure. Hence, long folds grad¬ 
ually become less probable and P{£/R) becomes more 
asymmetric due to relatively large populations of smaller 
folds. The evolution of P{£/R) with increasing R/r in 
Fig. IHKa) should similarly indicate the growth of spatial 
exclusion effects, which is counterintuitive since the pack¬ 
ing fraction (/> decreases with increasing the system size. 


The circular cross section of wires allows them to slide 
over each other rather easily (which results in slightly 
curved folds). Indeed, the strength of self-avoidance ef¬ 
fects is better captured by the length of the crumpled 
wire, rather than the volume excluded by it. Thus, we 
can ignore the role of wire thickness as a spatial exclu¬ 
sion constraint, and assign the entire contribution to the 
length of wire, reflected in the dimensionless quantity 
A=^~(A/r)i-5 (while ^=|^^(i?/r)-0-5). 

A closer look at the tail behavior of P{£/R) in Fig.|3l(b) 
verifies that the correlations grow with R/r. The experi¬ 
mental data is fitted to Gaussian, log-normal, and gamma 
distributions over the entire range of £/R. While the 
Gaussian function represents the distribution of random 
(uncorrelated) data, the gamma and log-normal distri¬ 
butions are associated with random events in the pres¬ 
ence of self-correlations or those that occur hierarchically 


25l-l28l| . The tail behavior is better captured by the 


Gaussian for small R/r, while there is a gradual crossover 
towards the two other functions as R/r increases. 

Self-avoiding random walk model — To better under¬ 
stand the mechanism of increasing spatial exclusion ef¬ 
fects, we model the folding process as a self-avoiding 
random walk (SAW) of the wire inside the confine¬ 
ment. While the existing SAW algorithms mainly follow 
stochastic Markovian dynamics to sample the ensemble of 
trajectories on regular lattices, here we propose an alter¬ 
native approach which accounts for the time evolution of 
the step size £. We first suppose that the strength of self¬ 
avoidance effects after the n-th fold is mainly controlled 
by the length of the inserted wire, i.e. the larger is the 
parameter A„=-^, the smaller is the success probability 
for the fold n-l-1 to be a long fold. The size £ of the next 
fold is obtained via the following algorithm: The distance 
£ between two randomly chosen points inside the sphere 
is taken as the trial fold size, and the proposal is accepted 
according to a Metropolis-like criterion with probability 


Pn+i{£)=Af ^ exp [-/3A„£/i?], 


(3) 


where A/’=^g^(l—exp[—/3A„£i„ax/£?]) is the normaliza¬ 
tion factor, and /? is a coefficient which depends on wire 
properties. In the case of rejection, a new £ is tried. Fi¬ 
nally, the cumulative length is updated as L„+i=L„-|-£ 
before starting the next step. Equation m assumes 
an exponentially lower acceptance chance for larger trial 
lengths £. Moreover, the acceptance probability decreases 
with increasing An, as the self-avoidance effects become 
more pronounced. The method indeed samples the fold- 
length landscape according to a Boltzmann-like distri¬ 
bution. Initially, the wire walks in free space (An=0), 
thus, Ri+i{£) = j^ independently of the trial length £. 
However, larger values of £ become gradually less prob¬ 
able with increasing A„. We perform extensive Monte 
Carlo simulations by adjusting £niax and the threshold 
value of (j) to the experimental data. The results shown 
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FIG. 4. (color online). The cumulative length Ln of the 
inserted wire after the n-th fold. The full (open) symbols 
indicate experimental (simulation) results. The simulation 
results are multiplied by 10 for clarity. Inset: Ln/R averaged 
over all experiments (solid line) and its standard deviation 
(shaded area). The dotted line is obtained from Eq. ([^. 


in Figs. [2J4] are in remarkable agreement with experi¬ 
ments (the straight folds in simulations lead to slightly 
lower values of {£)/R in Fig. [2{b) compared to curved 
folds in experiments). The exponents are not affected by 
the choice of /3, however, we treat it as a free param¬ 
eter to quantitatively reproduce the experimental data. 
We also checked that the power-law scalings can not be 
reproduced when replacing A„ with (/)„ in Eq. (|3]). Im¬ 
portantly, from the simulation results in Fig. EKb), we 
can specify the growing of self-correlations as dominant 
mechanism rather than hierarchical events. A numerical 
study of crumpled elastic papers also revealed that self¬ 
avoidance at high compression leads to self-correlations 
which modify the hierarchical folding process and the 
shape of the segment length distribution Q. 

The cumulative length of the inserted wire after the 
n-th fold qualitatively collapses onto a master curve for 
different values of R/r, as shown in Fig. ID The folds are 
initially independent of each other and Ln/R grows lin¬ 
early with n. The steps however become more correlated 
with increasing n, leading to a slower growth of Ln/R- 
From the scaling of L and N with R/r [Eqs. ([1]) and ([2|)] 
one obtains the asymptotic scaling L/i?^ 
which is consistent with experiments (see inset). Starting 
from (^1)=^^*^, the mean fold size at next steps can be 
estimated in terms of the accumulated wire length as 


^^max TD 

{£n+l)= Pn+l{i)idi= — 

Jo P'^n 


^_g/3An^niax/-R 


( 4 ) 


Hence, we obtain the following recursive analytical ex¬ 
pression for the injected length of wire after the n-th fold 


T f 0 \ 

^-max . X ' 


R 2R 


i=l 


R _ £ma.x/R 

J3L, l-exp[/3£inaxTi/i?2] 


( 5 ) 


in excellent agreement with the data [see Fig. lUinset)]. 

Buckling threshold — Finally, we clarify how the injec¬ 
tion of plastic wire eventually stops. According to the ex¬ 
perimental results shown in Fig. [SI the cutoff fold length 
ic (averaged over the last five folds) is independent of the 
system size R for a given wire radius. We suggest that 
the process stops when the fold size becomes so small 
that the applied feeding force F is unable to fold the 
wire anymore. The minimum threshold length can be 

obtained from the Euler buckling theory as £c=\J~^^'r‘^ 
for fixed-free end conditions, thus, tc is determined by 
the applied force and wire properties independent of R. 
For r=0.5 mm and a typical force of F=100A^, we obtain 
£c—6.0 mm, comparable to the experimental values. 

Figure |4] interestingly evidences a universal filling 
mechanism independent of the relative system size. The 
normalized cumulative wire length Ln/R versus the nor¬ 
malized length of the n-th fold £n/R collapses onto a uni¬ 
versal curve for different values of R/r (Fig. |S]). As the 
wire injection continues, £n/R gradually decreases until 
it reaches the minimum threshold value £c/R, where the 
process eventually stops. For a given value of r, is 
independent of A, thus, Ic/R decreases (i.e. shifts to the 
left in Fig. [5]) with increasing R. Consequently, the cutoff 
number of the steps at which the process stops increases 
with the system size. By calculating £c from the inser¬ 
tion force and wire properties, one can predict the total 
length of the crumpled wire for different system sizes. 

In conclusion, we studied the growing self-avoidance 
effects during the injection of plastic wires in frictional 
containers. The direct measurement of the evolution of 
fold-length distributions verified that the process goes 
on in a correlated way, such that the longer the injected 
wire, the stronger the spatial exclusion effects. The re¬ 
sults provide a new insight into underlying mechanisms 
of crumpling, which goes beyond a simple hierarchical de- 



FIG. 5. (color online). Gollapse of Ln/R vs. the scaled fold 
length £n/R for different values of i?/r in experiments (back¬ 
ground gray lines) and simulations (thick colored lines). Same 
colors as in Fig. (4) The symbols show the experimental cutoff 
values £c/R for r=0.5mm. Inset: £c versus the sphere radius 
R. The horizontal line shows the average value £c=5.8 mm. 
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scription of the process. The proposed sampling method 
of SAW trajectories would open the door to feasible 
Monte Carlo simulations of self-avoiding random walks 
in continuous space inside arbitrary confinements to ob¬ 
tain macroscopic quantities of interest, as an alternative 
to computationally expensive DEM methods. 

The experiments were done in the Complex Fluid Lab¬ 
oratory of lASBS. The authors gratefully acknowledge 
support from the lASBS Research Council and DFG 
within SFB 1027 (A7). 


* shaebani@lusi.uni-sb.de 
^ M.Habibi@uva.nl 

[1] P. K. Purohit, J. Kondev, and R. Phillips, Proc. Natl. 
Acad. Sci. U.S.A. 100, 3173 (2003); 1. Ali, D. Maren- 
duzzo, and J. M. Yeomans, Phys. Rev. Lett. 96, 208102 
(2006); E. Katzav, M. Adda-Bedia, and A. Boudaoud, 
Proc. Natl. Acad. Sci. U.S.A. 103, 18900 (2006). 

[2] J. Kindt, S. Tzlil, A. Ben-Shaul, and W. M. Gelbart, 
Proc. Natl. Acad. Sci. U.S.A. 98, 13671 (2001); D. E. 
Smith, S. J. Tans, S. B. Smith, S. Grimes, D. L. Ander¬ 
son, and G. Bustamante, Nature 413, 748 (2001). 

[3] J. H. Brackenbury, J. Zool. 232, 253 (1994). 

[4] H. Kobayashi, B. Kresling, and J. F. V. Vincent, Proc. 
R. Soc. B 265, 147 (1998). 

[5] G. A. Vliegenthart and G. Gompper, Nature Mater. 5, 
216 (2006). 

[6] T. Tallinen, J. A. Astrom, and J. Timonen, Nature 
Mater. 8, 25 (2009). 

[7] A. S. Balankin et ah, Phys. Rev. E 87, 052806 (2013); 
Physica A 388, 1780 (2009). 

[8] E. Sultan and A. Boudaoud, Phys. Rev. Lett. 96, 136103 


(2006). 

[9] Y.-C. Lin et ah, Phys. Rev. Lett. 103, 263902 (2009). 

[10] A. D. Gambou and N. Menon, Proc. Natl. Acad. Sci. 
U.S.A. 108, 14741 (2011). 

[11] H. Aharoni and E. Sharon, Nature Mater. 9, 993 (2010). 

[12] P. Lanzer, Mastering Endovascular Techniques: A 
Guide To Excellence (Lippincott Williams and Wilkins, 
Philadelphia, 2007). 

[13] R. Vetter, F. K. Wittel, and H. J. Herrmann, Nat. Com- 
mun. 5, 4437 (2014). 

[14] N. Stoop, J. Najafi, F. K. Wittel, M. Habibi, and H. J. 
Herrmann, Phys. Rev. Lett. 106, 214102 (2011). 

[15] J. Najafi, N. Stoop, F. K. Wittel, and M. Habibi, Phys. 
Rev. E 85, 061108 (2012). 

[16] M. Habibi, J. Najafi, and N. M. Ribe, Phys. Rev. E 84, 
016219 (2011). 

[17] J. A. Aguiar, M. A. F. Gomes, and A. S. Neto, J. Phys. 
A 24, L109 (1991). 

[18] P. G. de Gennes, Scaling Concepts in Polymer Physics 
(Gornell University Press, Ithaca, 1985). 

[19] M. A. F. Gomes, V. P. Brito, and M. S. Araujo, J. Braz. 
Ghem. Soc. 19, 293 (2008). 

[20] Y. Kantor, M. Kardar, and D. R. Nelson, Phys. Rev. 
Lett. 57, 791 (1986). 

[21] S. Havlin and D. Ben-Avraham, Phys. Rev. A 26, 1728 
(1982). 

[22] G. G. Donato, M. A. F. Gomes, and R. E. de Souza, 
Phys. Rev. E 66, 015102(R) (2002). 

[23] N. Stoop, F. K. Wittel, and H. J. Herrmann, Phys. Rev. 
Lett. 101, 094101 (2008). 

[24] F. Plouraboue and S. Roux, Physica A 227, 173 (1996). 

[25] A. J. Wood, Physica A 313, 83 (2002). 

[26] D. L. Blair and A. Kudrolli, Phys. Rev. Lett. 94, 166107 
(2005). 

[27] M. Adda-Bedia, A. Boudaoud, L. Bone, and S. Deboeuf, 
J. Stat. Mech.: Theory Exp. P11027 (2010). 

[28] S. Deboeuf, E. Katzav, A. Boudaoud, D. Bonn, and M. 
Adda-Bedia, Phys. Rev. Lett. 110, 104301 (2013). 



